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ABSTRACT 

We investigate properties of Galactic microlensing events in which a stellar object is lensed by a 
neutron star. For an all-sky photometric microlensing survey, we determine the number of lensing 
events caused by ^ 10® potentially-observable radio pulsars to be ~ 0.2 yr“^ for 10^° background 
stellar sources. We expect a few detectable events per year for the same number of background 
sources from an astrometric microlensing survey. We show that such a study could lead to precise 
measurements of radio pulsar masses. For instance, if a pulsar distance could be constrained through 
radio observations, then its mass would be determined with a precision of ^ 10%. We also investigate 
the time-scale distributions for neutron star events, finding that they are much shorter than had been 
previously thought. For photometric events towards the Galactic centre that last ~ 15 days, around 
7% will have a neutron star lens. This fraction drops rapidly for longer time-scales. Away from 
the bulge region we find that neutron stars will contribute ^ 40% of the events that last less than 
^ 10 days. These results are in contrast to earlier work which found that the maximum fraction of 
neutron star events would occur on time-scales of hundreds of days. 

Subject headings: gravitational lensing - pulsar: general - stars: neutron 


1. INTRODUCTION 

Neutron stars (NSs), as the densest observable objects 
in the Universe, allow us to study diverse phenomena in 
astrophysics and also probe fundamental physics at ex¬ 
tremely high density. For instance, determining whether 
pulsar-like compact objects are NSs or quark stars would 
have profound implications on the physics of condensed 
matter and on the nat ure of the strong interaction (e.g., 
Weber|[200^ p^[Mol for a review). 


Determining i\S masses is essential for constraining 
properties of their inner structure. The existence of 
high-mass NSs places strong constraints on the equation 
of state, while extremely low-mass NSs help distinguish 
self-bound quark stars from gravity-bound NSs. Recent 
precise high mass measurements of PSRs J0348-I-0432 
and J1614—2230 significantly constrained the equation of 
state of NSs and impr oved our understanding of pulsar- 


et al. 

2010| 

Ozel et al. 2010 

Lai&Xu|2011 

). However, up 

to now, all 

precise mass measurements are trom observa- 


tions of binary pulsars. Determining masses for isolated 
NSs is challenging, but knowing the mass distribution for 
such NSs is important as it may be different from binary 
NSs. 


to measure the masses of isolated NSs (e.g 

., fDai et al. 

2010 

Schwarz & Seidel||2002 

Horvath||1996 

I. There are 

two 

kinds of microlensing phenomena that 

occur when 


the line-of-sight to a distant object passes close to an 
intervening massive object. The first is known as photo- 
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metric microlensing, which is observed as the temporal 
brightening of the background star. For photometric mi¬ 
crolensing events the photometric magnification. A, de¬ 
pends on the angular separation between the lens and the 
source in units of the angular Einstein rad ius, u = 9 s\/9-e., 
as A ~ 1 -|- 2/it^ (Mao & Paczynski 1996). A typical an¬ 
gular Einstein radius tor a stellar rnicrolensing event in 
the Galaxy is approximately 1 mas. The second is known 
as astrometric microlensing and is observed as the angu¬ 
lar shift, S, of the apparent position of the background 
star. For angular separations larger than t he angular 


Einstein radius , 0 e, the shift S ~ 9b/ u (e.g., Belokurov 


& Evans 2002 hereafter BE02). Therefore, the cross- 


section for astrometric events will be much larger than 
for photometric events because the detectability of astro- 
metric events is not as strongly dependent on u. How¬ 
ever, it is currently impractical to carry out large area 
astrometric microlensing surveys. 

The Optical Gravitational Lensing Experiment 
(OGLE) 


(Udalski 


20031 is an ongoing photometric 
microlensing surveyl The third phase of this project, 
OGLE-III , monitored around 3 x 10® stars in the Galac¬ 
tic bulge (Szymahski et al.||2011 ), as well as fields in the 
disc and towards the IVIagellanic Clouds. The current 
phase of the OGLE project, OGLE-IV, is covering 
a significantly larger area and detecting many more 
events than OGLE-III. In the near future the Large 
Synoptic Survey Telescope (LSST) (for more detailed 
discussion see the LSST Science Book; |LSST Science| 
Collaboration et al.|2009 ) will also detect large numbers 


of microlensing events. 

There are currently no wide-area surveys for astromet¬ 
ric microlensing events and none planned in the fore- 
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seeable future. However, if an astrometric microlensing 
event could be predicted then a specihc observing cam¬ 
paign could monitor the particular area of sky in which 
the event will occur. It is possible to determine astro- 
metric parameters independently for NSs that can be 
observed as radio pulsars. For such pulsars it is there¬ 
fore possible to predict times of microlensing events. A 
further advantage is that the distance to the pulsar can 
be independently estimated. The properties of such a 
microlensing event would therefore be well constrained 
and would provide a direct measurement of the mass of 
the pulsar lens. Currently only a relatively small num¬ 
ber of radio pulsars are known. However, future ra¬ 
dio telescopes, such as the Fi ve-hundred-meter Aperture 

2000^ and the 


Spherical Telescope (FAST, Nan et al. 


Square Kilometre Array (SKA, Johnston et al. 2007), 


will discover a significant fraction ot all the potentially- 
observable radio pulsars in our galaxy and will be able to 
determine the astrometric parameters for many of these 
pulsars precisely. 

In the first part of this paper we study photometric NS 
leasing events. We calculate the expected event rates, 
the characteristic time-scale of NS events and determine 
the fract ion of total events caused by NSs. Previous stud¬ 
ies (e.g.,|Sartore & Treves|2010||Wood & Mao|2005l here¬ 
after WM05) were not based on the most up-to-date spa- 
tial and velocity distributions for the NS population (e.g., 
Hobbs et al.||2005[ Lorimer et al.||200^ |Faucher-Gigu^e| 
&: Kaspi||2006|). NS lensing events are only a small frac- 
tion ot the total number of events detectable in a lensing 
survey. We show that knowledge of the expected charac¬ 
teristic time-scale of NS events can signihcantly improve 
identification of such events. 

In the second part of the paper we make predictions 
for astrometric microlensing events caused by NSs. We 
provide expected events rates and determine the likely 
precision that could be achieved for NS masses. This has 
not been investigated previously (previous studies con¬ 
centrated on stellar lenses e.g., BE02). We conclude by 
discussing the implications of our results and providing 
predictions for future surveys. 

2. PROPERTIES OE PHOTOMETRIC 
MICROLENSING CAUSED BY NEUTRON 
STARS 

In this section, we investigate photometric microlens¬ 
ing events caused by NSs. NSs can lens emission from 
Galact ic stellar objects an d also from background galax¬ 
ies. As Tian & Mao (2012) have already studied NSs lens¬ 
ing faiiit background galaxies, we concentrate on Galactic 
background stellar objects. Our work follows the proce¬ 
dure described in WM05. In Section 2.1, we describe 
the spatial and velocity distribution of the NS popula¬ 
tion and the stellar models used in our simulations. We 
then present the event rates, time-scale distributions of 
microlensing events caused by NSs and their fractional 
contribution to the Galactic event rates. 

2.1. NS and stellar population distributions 

2.1.1. The distribution of the NS population 

We assume that the NS number density is proportional 
to that of radio pulsars. In cylindrical coordinates origi¬ 
nating at the Galactic Genter (GG), the NS number den¬ 


sity is expressed as 
p{R,z) = 


R^ 


i?0 + Rl 


exp 


-h 


R — -^0 

Rq + Rl 


X exp 


E 


(1) 


where i?i = 0.55 kpc, a = 1.64, b = 4.01, E = 330 pc 
and the distance from the Sun to the GG is set to be 
i?Q = 8.0 kpc. This assumes the radia l distri bution of 

(|20Mt 

-Cher-G: 


NSs as suggested by Yusifov fc Kiigtikl 
radial distribution was also used m [Faucher-Uiguere 


Such a 

1 


Kaspil (|2006D, and allows us to use their prediction for 
the number of potentially-observable pulsars (i.e. NSs 
that are radio-loud and beamed toward us). Eor the z 
distribution we use the scale height obtained by [Lorimer 


et al. 

TlTe 


(2006). 

tota. 


number of NSs in the Galaxy is still un¬ 
clear. Eollowing recent estimations (Keane & Kramer j 
2008t , we normalize the parameter A to give a total INS 
number of 10®. For potentially-observable pulsars A is 
normalized to give a total number of observa ble pulsars 
of 120,000 dFaucher-Giguere & Kaspi||2006|). We note 
that this number does not account fo r the population of 


Rotat ing Radio Transients (RRATs, McLaughlin et al. 


20061. Their inclusion wo uld significantly increas e the 


number of observable NSs (Keane &: Kramer||2008 ). 


We adopt a Gaussian model with a = 290 km s 


-1 

dFaucher-Giguere & Kaspi||2006|) for the distribution of 
the components ot the INS velocities. As will be defined 
in Section 2.2, the time-scale and event rate is inversely 
proportional and proportional to the lens-source rela¬ 
tive transverse velocity respectively, therefore increas¬ 
ing or decreasing the a of the NS velocity distribution 
will simply rescale the time-scale and event rate. Al¬ 
though an expon e ntial model is preferred in [Faucher^ 
Giguere & Kaspi (2006), the Gaussian model has also 
been shown to •g^l match the v elocity distribution of ra¬ 
dio pulsars ( Hobbs et al.||2005 l and provides a straight¬ 
forward distribution for us to use in our calculations. In 
the analysis below we use the Gaussian model by default, 
but for comparision we also give some results based on 
the exponential model. 

2.1.2. Stellar models 

In order to determine the properties of the background 
stars we consider stellar populations in the bulge and in 
the disk. For the bul ge distribution of stellar objects, we 
use the E2 model in iRattenbury et al. (20071. We t run¬ 


cate the bulge at a cor otation radius He = 3.5 kpc (Bis- 


santz & Gerhard 


2002). The model is normalized to the 


observed value ot the optical depth towards the Central 
Galactic R egion (Calchi Novati et al. 2008 Popowski 


et al. 2005). The disk distribution of stellar objects is 


relatively well co nstrained, and we u se the local vertical 


density model of Zheng et al. (|200I), as ext ended to the 


whole Galactic disk by 


We also need models 


(2003). 


"Han & Gould 

ot the lensing objects (both NSs 
and other objects). For the lens mass function, we follow 
an identical procedure to that described in WM05. This 
is based on a two part power-law distribution. 


dN _ ^ / M y 
dM “ ^ I MbrkJ 


(2) 
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where Mbrk = 0.7 Mq, a = —2.0 for M > Mbrk, and 
a = —1.3 for M < Mbrk- 

We adopt the kinematic model of WM05, which de¬ 
scribes the velocities of the lenses, sources and the ob¬ 
server. The lens-source relative transverse velocity, v, 
is calculated according to Eq. 7 and 8 in WM05. The 
conversion of velocities from the Galactocentric cylindri¬ 
cal coordinates int o solar-centric spherica l coordinates is 
given by Eq. 15 in Wang & Smith (2011). 


To calculate the event rates and mean time-scales, we 
use the adaptive Monte Carlo integration algorithms in¬ 
corporated in the GNU Scientific Library^. In order to 
improve the calculation precision we use a sampling num¬ 
ber of 10^ over the multi-dimensional volume when we 
carry out the Monte Carlo integrations. 

2.3. Photometric microlensing results 
2.3.1. Event rates 


2.2. The photometric microlensing model 

The time-scale of a microlensing event is the time taken 
for a source to cr oss the Einstein radius of the lens, te, 
and is defined as (Paczynski||1996|) 


= 


rE 


rE = 


UGM Jd(-Ps --Pd) 
c2 ^ D. 


( 3 ) 


where G is the gravitational constant, c is the speed of 
light, Ds is the source distance and is the lens dis¬ 
tance. 

The event rate, T, is defined as the number of photo¬ 
metric microlensing events per unit time for a given num¬ 
ber of background sources. The event rate depends upon 
the background source number density, ps{Ds), and the 
mass density of lensing o bjects, Pm(-Od). Following Ki- 
raga & Paczynski ( 1994|) and WM05, we assume that 
the traction ot stars brighter than some luminosity L is 
proportional to and then the event rate is expressed 
as: 


P =- 


4G1/2 


c 

rD; 




/o ° dDdPm(Dd)t^[Dd(Ds - Dd)/MD,]i/^ 

/“dD,D2+2^Ps(Ds) 


( 4 ) 


where v is the lens-source relative transverse velocity. 

We use P to model numerous phenomena that are not 
explicitly included in the equation. To account for the 
luminosity function of the background sources, the op¬ 
tical survey sensiti vity and the Galactic exti nction, we 
set —3 < P < —1 (Kiraga & Paczynski 1994). We also 
calculate our results with /I = U, which corresponds to 
the maximum possible event rate. This is a simplified 
approach, but has been applied by previous studies on 
the Galactic microlensing event rate and is valid as we 
only report basic properties of NS lensing events and 
do not attempt to predict the exact number of events 
detected with a particular set of survey parameters. Al¬ 
though extinction will clearly affect our estimates for the 
event rates, implementing this in a consi stent manner 
is a v ery complex task (see, for example, Kerins et al. 
2009). However, since extinction results iii a change in 
the slope of the luminosity function, this can be reflected 
in the choice of the P parameter, i.e., large negative val¬ 
ues of beta can be used to account for high extinction. 
We consider stellar objects in the Galactic disk and in 
the bulge as background sources. As the disk stars can 
have small distances, setting —3 < P < —1 will result in 
unreasonably small event rates. We account for this by 
always setting /3 = 0 for disk sources within 4.5 kpc from 
the observer. 


Our results give a maximum possible all-sky averaged 
event rate (P = 0) of 4.2 x 10“^ yr“^ for NS lenses. 
Towards the GG and the Baade’s window (BW)^ re¬ 
gion this maximum event rate is 6.7 x 10“^ yr~^ and 
3.6 X 10“”^ yr“^, respectively. Setting P = —1,—2,—3, 
the all-sky averaged event rate is 1.3 x 10“^ yr”^, 
7.5 X 10“® yr“^ and 7.4 x 10“® yr“^, respectively. In Ta¬ 
ble!^ we present the all-sky averaged event rate, rates for 
I 6 f< 5° and | 5 | > 5° and rates towards the GG and BW, 
for P = 0,—1,-2,—3 for 120,000 potentially-observable 
radio pulsars. Since we have assumed that the popula¬ 
tion model for radio pulsars is identical to that of NSs, 
except for the normalization, this means that the event 
rates for potentially-observable radio pulsars is a factor of 
10®/ (l.2 X 10®) smaller than that of NSs. Using the ex- 
ponential model proposed by |Faucher-Giguere & Kaspij 
|200^ , we obtained an all-sky averaged event rate for 
potentially-observable radio pulsars of 4.2 x 10“^® yr“^ 
{P = 0), compared with 5.Ox 10“^^ yr“^ for the Gaussian 
model. 

We show in Fig. the expected NS event rates in the 
central Galactic region. The contour levels represent the 
event rate using a logarithmic scaling, log(r/10“® yr“^). 
Two sets of contours are shown. Contours for P = —1 are 
shown with solid lines. Contours for P = —2 are shown 
as shaded regions. As expected, the event rate is higher 
on the Galactic plane and decreases at higher Galactic 
latitudes. As P varies from —1 to —2, the distribution 
of the event rates becomes flatter, and the rate decreases 
much more significantly on the Galactic plane compared 
to that at a higher Galactic latitude region. 

The reduction in event rates as P decreases and the 
flattening of the event rate distribution can be explained 
as: 

• Decreasing P implies that fewer distant sources 
contribute to the event rate. Therefore the total 
event rate will decrease as P decreases. 

• Towards the GC the stellar density is higher than 
in the Galactic disk. Decreasing P therefore has a 
larger effect in the GC region than elsewhere. 

• Sources in the bulge have a high velocity disper¬ 
sion. Therefore decreasing P effectively removes 
more high-velocity stellar objects in the bulge com¬ 
pared with the disk. 

• Sources on the far side of the Galactic disk rotate 
in the opposite direction to the observer. These 
sources therefore contribute a significant number of 


^ http://www.gnu.org/software/gsl/ 

^ (L, b) = —2.75“), L and b are the Galactic longitude and 

latitude. 
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Fig. 1.— Map of the microlensing event rates caused by NSs, 
towards the central Galactic region. Contours with solid lines rep¬ 
resent results for = — 1. Contours filled with colors represent 
results for j3 = —2. OGLE-III fields are shown with diamond sym¬ 
bols. Contour levels represent the event rate in logarithmic scaling, 

log(r/10“^ yr“^)- 


events as their relative velocities to the NS lenses 
are high. Decreasing j5 removes these distant back¬ 
ground sources. 


OGLE-III fields^ are shown in Fig. wit h diamond 


symbols. F or 340 x 10® background sources (Szymahski 


et ari[2011|) and taking /3 = —1, we predict there should 
To 


have been around 50 NS events per year and about a 6% 
chance that there has been an event from a potentially- 
observable radio pulsar during the 9-year duration of the 
survey. The current phase of the OGLE project (OGLE- 
IV) started in 2010 and this is now detecting around 
four times as many microlensing events per year in the 
bulge compared to OGLE-IlC. This could therefore de¬ 
tect around two hundred NS events per year and, given 
the event rate for potentially-observable radio pulsars, 
there is a ~ 23 per cent chance it will detect one by the 
end of the decade. For future deeper and wider optical 
surveys the prospects for discovering pulsar microlensing 
events are promising. For example, LSST is expec ted to 
monitor ^ 10^® background Milky Way sources (Ivezic 


et al. 


20121. Taking /3 = — 1, we predict there would be 


lUUl) events caused by NSs per year and ~ 1.5 events 
caused by potentially-observable pulsars per decade. 

Note that our photometric event rates are likely to un¬ 
derestimate the actual number of observable events. This 
is because our calculations are based on the definition of 
a microlensing event as one for which the lens passes 
within one Einstein radius of the source, corresponding 
to an amplification of 0.32 mag. However, in practice it 
is possible to detect events with much weaker amplifica¬ 
tion. If we instead define an event as one for which the 
amplification is 0.1 mag, then the event rate will be ~ 70 
per cent higher. There is also the fact that the actual 
number of source stars may be underesti mated due to 


blend ing from unresolved source stars (e.g. Smith et al. 


20071. 


^ data from http://oKle.astrouw.6du.pl/ 
http://ogle.astrouw.edu.pl/ogle4/ews/ews.html 


TABLE 1 

The event rates caused by potentially-observable radio 
PULSARS, IN UNITS OF 10“^^ yr“^- TO OBTAIN THE EVENT RATE 
FOR NS LENSES, THE NUMBERS IN THIS TABLE SHOULD BE 
MULTIPLIED BY A FACTOR OF 10®/ (l.2 X 10®). 
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Fig. 2. — Microlensing event rates caused by NSs as a function 
of time-scale. The solid, dotted and dashed lines represent the all¬ 
sky averaged event rate and the event rates towards the GC, BW, 
respectively. In comparison, the time-scale distribution for NSs 
towards the BW based on WM05’s model is shown as indicated. 



Fig. 3. — Microlensing event rates caused by NSs as a function 
of time-scale towards the GC. The solid, dashed, dotted and dot- 
dashed lines represent event rates for /3 = 0,-1,—2,—3, respec¬ 
tively. 


2.3.2. Time-scale distributions 

By restricting the Monte Carlo integrations in a pa¬ 
rameter space corresponding to a certain time-scale, we 
can directly obtain the event rate caused by NSs as a 
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function of time-scale. These are shown in Fig. for 
/3 = 0. The time-scale distributions for the all-sky av¬ 
eraged event rate is shown as the solid line. The event 
rates towards the GC and BW are presented as dotted 
and dashed lines respectively. We note that the distribu¬ 
tions have asymptotic power-law tails fol l owing the ex¬ 
pected form given by Mao & Paczynski (19961. In all 
cases the event rate peaks at a time-scale ot 12 days. 
The mean time-scales for the all-sky averaged, GG and 
BW event rates are ^ 19.7, 19.3 and 19.4 days respec¬ 
tively. For comparison we also overlay the prediction 
from WM05 towards the BW which has a mean time- 
scale of 47.3 days. Our new results therefore imply that 
the majority of NS events will have shorter time-scales 
than predicted by WM05, which is mainly because of 
the relatively high velocity of NSs. Using the exponen¬ 
tial model, we obtained a mean time-scale for the all-sky 
averaged event rate of ^ 23.7 days. 

In Fig. 1^ we show the effects of changing (3 on events 
towards tli^e GC. The difference in the event rate be¬ 
tween /3 = — 1 and —2 is much larger than for other /3 
values because the background source density starts to 
decrease with increasing distance according to the 
term in Eq. ffl Below the mean time-scale the choice of 
P has little effect on the results. /3 does affect the longer 
time-scales with, as expected, smaller event rates as /3 
decreases. However, the mean time-scale is slightly de¬ 
pendent on P with the mean time-scales being 19.3, 17.9, 

13.5 and 13.2 days for /3 = 0 to —3 respectively. This 
shift was not initially expected as decreasing P reduces 
the averaged transverse velocity between the sources and 
lenses leading to an expected decrease in the number of 
short-duration events. In order to understand this result, 
we calculated the Einstein radii averaged over all lenses 
and the averaged transverse velocity between the lenses 
and sources. We found that, as p decreased, the averaged 
Einstein radii decreases faster than the reduction in aver¬ 
aged transverse velocity. This implies that the resulting 
lensing events will occur by lenses with smaller Einstein 
radii and therefore the time-scale of these events will be 
shorter. 

WM05 predicted the averaged time-scales of all (not 
just NS) Galactic events as a function of Galactic co¬ 
ordinates. The short duration events that we are pre¬ 
dicting for NS lenses will lower the mean time-scale of 
these events. Our map is presented in Fig. and can be 
compared with Fig. 3 in WM05. The solid contours rep¬ 
resent P = —1 and the shaded regions as P = —2. The 
OGLE-III fields are shown with diamond symbols. Our 
averaged time-scales towards the BW are 23.3, 23.9 and 

36.5 days for /? = 0, — 1, —2, respectively. For P = 0,-1, 
our results are similar to, but slightly smaller than, the 
corresponding results in WM05 and generally agrees with 
the a ctual measurement of 28.1 ± 4.3 days (Sumi et al. 
20061 from the OGLE survey. 


2.3.3. Fractional contribution to the Galactic event rate 

We have shown that NS ev ents have relatively short 
time-scales. Sumi et al. (20111 showed that Jupiter-mass 
objects have even shorter tiiiie-scales, while other stel¬ 
lar objects can have longer time-scales. In this subsec¬ 
tion, we consider the fraction of Galactic events that are 
caused by NS lenses for different time-scales. 

In Fig. we present the fraction of events that are 



10 5 0 -5 -10 

Galactic longitude (deg) 


Fig. 4.— Maps of the averaged time-scales of Galactic microlens¬ 
ing events. Contours with solid lines represent results for /3 = —1. 
Contours filled with colors represent results for (3 = —2. Contour 
levels represent averaged time-scales in unit of days. OGLE-III 
fields are shown with diamond symbols. 


caused by NSs as a function of time-scale. The solid 
line refers to events towards the BW and the dashed line 
towards the GC. We also show results for two off-GC di¬ 
rections, {l,b) = (—20°,0°) and {l,b) = (0°,—10°). For 
(/, b) = (—20°, 0°) (on the disk, but away from the bulge) 
NS lenses contribute up to ^40% of the events that last 
less than ~ 10 days, but we note that the total number of 
events in this direction is low. For the three other sky po¬ 
sitions, the fractional contribution exhibits a clear peak 
around a time-scale of ~ 15 days and decreases at higher 
and lower time-scales. Our results do not agree with ear¬ 
lier predictions from WM05 (shown as a dot-dashed line 
in the Figure). Their work led to very small fractional 
contributions from NSs at ~15 days and much larger frac¬ 
tional contributions on time-scales of months to years. 
On the contrary, in our calculations we have found that 
NSs mainly contribute to short-duration Galactic events 
because of their relatively high velocities (a factor not 
considered by WM05). Away from the Galactic bulge re¬ 
gion, as the rate of short-duration Galactic events drops, 
the fractional contribution from NSs increases rapidly at 
short time-scales. 

In Fig. we show the fractional contribution towards 
the GC from NSs for different values of p. In all cases 
the fractional contribution peaks at ~ 10 days and we 
cannot reproduce the WM05 results by modifying p. As 
P decreases the fractional contribution on time-scales ^ 
10 days increases. For P = —3 almost all events at this 
time-scale would be caused by NS lenses. 

The fractional contribution from NS lenses changes as 
a function of sky position. Fig. shows the fractional 
contributions from NSs in the central Galactic region. 
Contours indicated with solid lines are for /3 = — 1. Con¬ 
tours shown as shaded regions are for p = —2. The 
OGLE-III fields are indicated using diamond symbols. 
The fractional contribution from NSs clearly varies with 
Galactic latitude with higher contributions from NSs oc¬ 
curring away from Galactic bulge. 

Averaged over all time-scales and over the entire sky 
we find that NSs contribute ~ 5.4% {P = 0) of the total 
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Fig. 5.— Fractional contributions to the total Galactic event rate, 
as a function of the event time-scale, from NSs. The solid, dashed, 
dotted lines and dot-dashed line show the contribution from NSs 
towards the BW, GC and off-GG directions {l,b) = (— 20°,0°) 
and (l,b) = (0°,—10°), respectively. In comparison, the fractional 
contribution from NSs towards the BW based on WM05’s model 
is shown as indicated. 



Fig. 6. — Fractional contributions to the total Galactic event 
rate, as a function of the event time-scale, from NSs, towards the 
GG. The solid, dashed, dotted and dash-dotted lines represent the 
contributions for /3 = 0, —1, —2, —3, respectively. 


number of Galactic events. In Table 2, we present the 
fractional contribution to the total averaged event rate 
from different types of lensing object. The different ob¬ 
jects considered are brown dwarfs (BDs), main-sequence 
stars (MSs), white dwarfs (WDs) and black holes (BHs). 
We also show contributions for different values of (3. For 
/3 = 0, the NS contribution overall is small, however for 
other values of /3 the contribution from NSs is much more 
signihcant. By increasing and decreasing the a of the NS 
velocity distribution by 30%, we got a fractional contri¬ 
bution from NSs of approximately 6.5% and 3.7% re¬ 
spectively (/3 = 0). Using the exponential model, the 
fractional contribution from NSs is approximately 4.2% 
(^ = 0). 


TABLE 2 

Percentage contributions to the total event rate, from 
DIFFERENT TYPES OF LENS. RESULTS FOR /3 = —1, —2, —3 ARE 
also PRESENTED. In COMPARISON, RESULTS BASED ON WM05’S 
MODEL ARE ALSO SHOWN. 


0 

BD 

Types of lens 
MS WD BH 
(per cent) 

NS 

0 

16.7 

60.5 

16.5 

0.9 

5.4 

-1 

15.8 

57.2 

15.5 

0.9 

10.6 

-2 

10.6 

38.2 

10.4 

0.5 

40.3 

-3 

10.0 

36.0 

9.8 

0.5 

43.7 

0 (WM05) 

17.2 

62.1 

16.9 

0.9 

2.9 



Galactic longitude (deg) 


Fig. 7. — Maps of the percentage contribution to the total Galac¬ 
tic microlensing event rates from NSs. Gontours with solid lines 
represent results for /3 = — 1. Gontours filled with colors represent 
results for (3 = —2. Contour levels represent the percentage con¬ 
tribution in per cent. OGLE-III fields are shown with diamond 
symbols. 


3. MASS MEASUREMENT THROUGH 
ASTROMETRIC MICROLENSING 

We have investigated the properties of photometric 
microlensing events caused by NSs. However, the pho¬ 
tometric microlensing event rate caused by potentially- 
observable radio pulsars is only ^ 2 per decade for 10^° 
background sources. As astrometric microlensing has a 
larger cross-section than photometric microlensing, more 
astrometric microlensing events are expected to be de¬ 
tected. It is not currently feasible to carry out large scale 
astrometric microlensing surveys. A more promising way 
to search for events caused by observable radio pulsars is 
to use the known pulsar astrometric parameters to make 
an initial guess for when and where microlensing events 
will occur. This approach can also be applied to photo¬ 
metric microlensing, but the larger cross-section means 
that it is especially promising for astrometric detections. 

The method we used to estimate the photometric mi¬ 
crolensing event rate is also valid for astrometric lens¬ 
ing, but the larger cross-section must be accounted for. 
The cross-section and the amplitude of an astrometric 
microlensing event can be characterized by the impact 
parameter, uq, which is defined as the smallest angular 
separation between the lens and the source in units of 
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the angular Einstein radius. The astrometric microlens¬ 
ing event rate will simply be our determination of the 
photometric event rate scaled by the increase in the cross- 
section. For instance, if the cross-section is ^ 10 times 
larger (uq = 10, discussed in more detail in Section 3.2) 
then, for 120,000 potentially-observable radio pulsars, 
the all-sky averaged event rate with /3 = — 1 will then 
be ^ 1.5 per year for 10^° background sources. 


3.1. Astrometric microlensing model 

BE02 presented methods to determine the precision 
that could be achieved for stellar lens mass estimates ob¬ 
tained through astrometric microlensing. They showed 
that there was a total of 11 parameters to completely 
parameterize such events. We apply their methods to 
determine the precision achievable for observable radio 
pulsars. 

For radio pulsar lenses, the positions, proper motions 
and distances of lenses can, in principle, be measured or 
estimated through radio observations. Pulsar positions 
and proper motions can be determined either through 
radio pulsar timing or interferometric observations. For 
millisecond pulsars that have been observed over many 
years, the precision of positions and proper motions de¬ 
termined through pulsar timing c ould reach < 1 rnas an d 
< Imasyr”^, respectively (e.g., Verbiest et al. 2008). 


For pulsars whose positions and proper motions cannot 
be obtained precisely through pulsar timing, interfero¬ 
metric observations can be used^. In our model, we con¬ 
sider an ideal case in which we assume that the radio 
pulsar lens positions and proper motions are known pa¬ 
rameters and we neglect their uncertainties. The impli¬ 
cations of these assumptions are discussed at the end of 
Section 3.2. 

Pulsar distances are much harder to measure than po¬ 
sitions or proper motions. For some pulsars precise dis¬ 
tances can be obtained from measurements of parallax. 
For most pulsars distances can only be estimated from 
dispersion measures. The conversion from dispersion 
measure to distance depends on a models for the dis¬ 
tribution of free electrons in the Galaxy. As such models 
are not very well constrain ed this results in distanc e un¬ 
certainties of around 20% (|Taylor & Cordes||1993|). Al¬ 
though models of the Galactic tr ee electron density ar e 
continually being improved (e.g., Gordes fc Lazio|[2002 ), 
it is not clear whether distances will be measured suf- 
ficiently for most pulsars at the time that astrometric 
lensing events from pulsars will be detected. Therefore, 
in our analysis we discuss two models. In the first model, 
which we label as Model Up"™", we assume that the 
pulsar distance is known. In the second model. Model 
Hp",.*^"™" unknown we assume that it is not. 

For Model six independent parameters need 

to be computed from the data. These are the position 
(0s,o), proper motion (/Xg), and parallax (tTs) of the back¬ 
ground source and the mass (M) of the lens. We note 
that it is possible that the astrometric parameters for the 
background source is also precisely known. If true, then 

® We note that the dynamical solar system ephemeris used for 
pulsar timing is different from the Inte rnational Celestial Reference 
Frame (ICRF). [Madison et al.|(|2013f computed a transformation 
between the ICKh' and pulsar timing frames and quantitatively 
discussed how the transformation will improve in the coming years. 


the only unknown parameter will be the mass. For this 
work, we assume that these parameters do need to be 
determined from the data. For Model 13™,,^"™", we have 
an additional parameter, the pulsar parallax (tti). 

As an example, we consider a specific pulsar, 
PSR B1929-I-10. This pulsar has a large proper mo¬ 


tion (103 mas yr Ghatterjee et al. 2004) and lies in 


close projection to a number of potential background 
source stars, with seven objects within 5 arcsecs accord¬ 
ing to our catalog comparison between known radio pul¬ 
sars® ( Manchester et al.|2005 1 and source tables from the 


UKIRT In frared Ueep Sky Survey (UKIDSS; La wrence 
et al.|200'7|). The distance to this pulsar is 0.36 kpc (|Gh^ 


terjee et al. 2004). We assume that this pulsar is likely 


to pass one of the background sources and could be stud¬ 
ied with an astrometric monitoring campaign that carries 
out observations of the background source with a uniform 
observing cadence over a period of four years (centered 
on the time of closest angular separation between the 
lens and the source). Following BE02, we assume that 
the background source is at a distance of 5 kpc with a 
proper motion of 7 mas yr“^. 

We note that our selection of lens and source proper 
motions and distances will not limit the generality of 
our analysis, since the amplitude of an astrometric mi¬ 
crolensing event is only characterized by Uq. However, 
the uncertainty on the mass measurements will depend 
on various parameters that need to be assumed in our 
modeling. The first is the pulsar mass. In our analy¬ 
sis, we use M = 1.4 Mq as a typical pulsar mass, but 
also trial different values. We also require an estimate 
of the astrometric accuracy that can be achieved for the 
background stellar object. We trial two possible values, 
cr = 150 and 450/ras. a = 150/ras corresponds to a rel¬ 
atively bright star at around V = 16 mag for the Gaia 
satellite^. Future telescopes such as the Wide-Field In¬ 
frared Survey Telescope (WFIRST) have the potential to 
obtain astrometric prec ision of a few tenth s of a mas for 
even fainter stars (e.g., Spergel et al.||2013 |. Finally, the 
uncertainty also depends on the nurnber of epochs that 
are observed during the four year monitoring program. 
In the results below we trial a range of possible values. 


3.2. Mass determination results 

Fig. [8] shows the uncertainties in M as a function of 
the dimensionless impact parameter uq for the two mod¬ 
els (in which the pulsar distance is known or not). We 
present results for our two models with different sam¬ 
plings (for 15 and 100 observation epochs). There is 


® http://www.atnf.csiro.au/research/pulsar/psrcat 
^ The end-of-mission parallax error for such a star is around 
20-40 fi&s. The single epoch positional precision is significantly 
worse. We estimate how much worse using information provided 
on http://www.cosmos.esa.int/web/gaia/science-performeince 
as follows, 'i'he positional precision for a single epoch measurement 
is more precise than the parallax precision by a factor of 0.743, but 
the single-epoch precision is a factor of 4.3 times worse than that 
available at the end-of-mission. From this we estimate that the 
single epoch positional precision is around three times worse. We 
also need to incorporate a factor of y/2 to account for the fact at 
each epoch Gaia only measures the one-dimensional position along 
the scanning direction. This gives an effective two-dimensional as¬ 
trometric precision of around 90 to 180 //as. Note that for faint 
stars {V = 20 mag) the end-of-mission parallax precision drops to 
around 130 to 600 /ias, meaning that the effective two-dimensional 
astrometric accuracy will be around 0.5 to 3 mas. 
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Fig. 8.— The percentage error in estimation of M as a function 
of impact parameter uq. We set M = IAMq and cr = 150^as. 
The solid and dashed lines represent the errors with 15 and 100 
observation epochs, respectively, and pulsar distances are supposed 
to be measured precisely (Model While the dotted and 

dot-dashed lines represent the errors for which pulsar distances are 
set to computed from the data (Model 



Fig. 9.— The percentage error in estimation of M as a function 
of the number of observation epochs, for Model Results 

for different uq and cr are shown. 


a large variation in the ability to determine the pulsar 
mass. The uncertainties in M increase with increas¬ 
ing Uq and decrease with the number of epochs. If the 
pulsar distance is known then its mass could be deter¬ 
mined at least two times better than if the mass is not 
known. Note that for Gaia the average number of epochs 
is around 70 (|de Bruiine||2012|). 

Fig. [9] and b ig. |l(j| show the uncertainties in M as a 
function of the number of observation epochs. In Fig. 
we assume that the pulsar distance is known, while in 
Fig. we assume it is unknown. For both figures, 
we present results for uq = 10 and uq = 30, and for 
tr = 150 fias and a = 450 ^as. For more than ~ 100 
epochs the uncertainty is approximately constant. More 
frequent sampling is therefore not necessary and will not 
improve the mass determination. As the astrometric pre¬ 
cision improves from 450/xas to 150 /ras, the uncertainty 
in M decreases by more than a factor of two. Clearly, it 
is therefore essential that the background source is mon¬ 
itored wit h h igh precision astrometry. 

In Fig. 0 we set different pulsar masses and show 
the uncertainties in M as a function of the number of 
epochs for Model In addition to a typical pulsar 

mass of 1.4 Mq, we also show results for an extremely 
low-mass case, 0.4 Mq, and a relatively high-mass case, 
2.4 Mq. We set ug = 10 and a = 450//as. As the pulsar 
mass increases the uncertainties become smaller. For the 
high-mass case, M = 2.4 Mq, the uncertainty is six times 
smaller than that of the low-mass case, M = 0.4 Mq. 

In the analysis above, we have presented results for an 
ideal case in which we treat pulsar positions and proper 
motions as precisely determined values and we neglect 
their errors. If the pulsar astrometric parameters are 
poorly known then those parameters would need to be 
included in the model when determining the lensing event 
properties. This would significantly degrade the mass de¬ 
termination precision. Here, we investigate an intermedi¬ 
ate case in which the astrometric parameters have been 



Fig. 10.— The percentage error in estimation of M as a function 
of the number of observation epochs, for Model . Results 

for different uq and cr are shown. 


measured, but do have small uncertainties. We there¬ 
fore model them as additional astrometric uncertainties 
of background sources and then apply the same method 
described above to estimate the uncertainty in pulsar 
masses. We assume that the uncertainties on the pulsar 
position and proper motion measurements are 0.01 mas 
and O.lmasyr”^ respectively. This is similar to that 
obtainable with current Very Long Baseline Interferom¬ 
etry (VLBI) results (e.g., Deller et al.| 20091. For Model 
assuming ug = 10 and 50 epochs, the uncer¬ 
tainty on the mass measurement increases slightly from 
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Fig. 11.— The percentage error in estimation of M as a function 
of the number of observation epochs, for Model Results 

for the pulsar mass of 0.4, 1.4 and 2.4 Mq are shown. 

10% to 12% compared with our earlier analysis. For 
larger pulsar position and proper motion errors of 0.1 mas 
and lmasyr“^, the uncertainty of mass measurement 
increases to 22%. Therefore, we conclude that the er¬ 
rors of pulsar position and proper motion can introduce 
additional uncertainties to pulsar mass measurements, 
but the effect is relatively small and can be significantly 
reduced by precise pulsar position and proper motion 
measurements already achievable with current VLBl ob¬ 
servations. 


4. CONCLUSIONS 


A number of current and future surveys are carrying 
out photometric monitoring campaigns of the Galactic 
bulge and d isc, such as OGLE-I V, Vista Variabl e in the 

^1 ' 


Via Lactea (Minniti et al.|2010) and WFI RST (Spergel 


et al. 20131. In the near future the LSST (Ivezic et al. 
2UUSp will monitor ~ 10,000 square degrees of sky us- 


mg pairs of 15-second exposures twice per night every 
three nights on average. This will lead to a typical 5cr 
depth for point sources of r ~ 24.5 (AB). Currently it 
is unclear to what extent LSST will survey the Galactic 
disc and bulge. Even without large portions of the bulge 
and disc, LSST will monitor 10 billion Milky Way stars 
and, if a strategy is chosen to include the disc and bulge, 
then this number will increase significantly. A specific 
3-year HST program (192 orbits) to detect microlensing 
events caused by non-luminous isolated black holes and 


proposed ( 

Sahu et al. 

2012) 

Uur modeling has s 

nown 


• The time-scale of microlensing events caused by 
NSs is much shorter (~ 20 d) than previously 
thought. 


• Events near the Galactic centre that have a time- 
scale of ~ 15 d will contribute to approximately 7% 
of the total number of events. Away from the bulge 
region the fractional contribution, for events with 
time-scales of < 20 d will increase to around 40%. 


• For an astrometric microlensing event caused by a 
known radio pulsar, pulsar masses could be deter¬ 
mined with a precision of ~ 10% if pulsar distances 
are known. With a sufficient number of observa¬ 
tions of the event, the mass could be measured to 
^ 25% if the distance is not known. 


For ongoing and upcoming surveys, we recommend: 

• To maximize NS events in the microlensing dataset, 
surveys of the Galactic disc and bulge should be 
carried out. 


• NS events should be identified by searching for 
short duration (~ 20 d) events (instead of long du¬ 
ration events that had been previously considered) 
and radio telescopes should observe the possible 
candidates in order to identify radio pulsars. 

• As new pulsars are discovered, their astrometric pa¬ 
rameters should be determined precisely and com¬ 
pared with catalogs of stellar positions. If a leasing 
event is likely then an astrometric leasing moni¬ 
toring campaign should be carried out around the 
predicted time of the event. 


• Although these targeted campaigns are clearly the 
best approach to detect astrometric microlensing 
events from potentially-observable ra dio pulsars. 


the p rospects for ESA’s Gaia mission (de Bruijne 


2012) are not negligible. In total Gaia will observe 


T(r stars and so, according to our estimates, this 
equates to around one such event over the lifetime 
of the mission (five years). 


In the more distant future, large radio telescopes will 
have discovered most of the potentially-observable ra¬ 
dio pulsars in the Milky Way. It is therefore likely that 
the mass of isolated NSs will be measurable through mi¬ 
crolensing events leading to improved knowledge of the 
pulsar population. Extremely high-mass and low-mass 
pulsars can then be studied in detail to provide new in¬ 
sight into the NS interior. 
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